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Fiber-reinforced  polymer  matrix  composite  materials  display  quite  complex  deformation  and  failure 
behavior  under  ballistic/blast  impact  loading  conditions.  This  complexity  is  generally  attributed  to  a 
number  of  factors  such  as  (a)  hierarchical/multi-length  scale  architecture  of  the  material  microstructure;  (b) 
nonlinear,  rate-dependent  and  often  pressure-sensitive  mechanical  response;  and  (c)  the  interplay  of  various 
intrinsic  phenomena  and  processes  such  as  fiber  twisting,  interfiber  friction/sliding,  etc.  Material  models 
currently  employed  in  the  computational  engineering  analyses  of  ballistic/blast  impact  protective  structures 
made  of  this  type  of  material  do  not  generally  include  many  of  the  aforementioned  aspects  of  the  material 
dynamic  behavior.  Consequently,  discrepancies  are  often  observed  between  computational  predictions  and 
their  experimental  counterparts.  To  address  this  problem,  the  results  of  an  extensive  set  of  molecular-level 
computational  analyses  regarding  the  role  of  various  microstructural/morphological  defects  on  the  Kevlar® 
fiber  mechanical  properties  are  used  to  upgrade  one  of  the  existing  continuum-level  material  models  for 
fiber-reinforced  composites.  The  results  obtained  show  that  the  response  of  the  material  is  significantly 
affected  as  a  result  of  the  incorporation  of  microstructural  effects  both  under  quasi-static  simple  mechanical 
testing  condition  and  under  dynamic  ballistic-impact  conditions. 
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1.  Introduction 

The  subject  of  the  present  study  is  high  specific-strength  and 
high  specific-stiffness  polymeric  fiber-reinforced  (e.g.,  Kevlar®, 
Twaron®,  etc.)  polymer-matrix  composite  materials  and  struc¬ 
tures.  These  materials/structures  are  commonly  used  in  various 
protective  systems  whose  main  requirement  is  a  high  level  of 
mass-efficient  penetration  resistance  against  high  kinetic  energy 
projectiles  (e.g.,  bullets,  mine,  IED  or  turbine  fragments,  etc.). 
Development  of  the  aforementioned  protection  systems  is 
traditionally  carried  out  using  legacy  knowledge  and  extensive 
fabricate-and-test  procedures.  Since  this  approach  is  not  only 
economically  unattractive,  but  is  often  associated  with 
significantly  longer  lead  times,  it  has  gradually  become 
complemented  by  the  appropriate  cost-  and  time-efficient 
computer-aided  engineering  analyses.  This  trend  has  been 
accelerated  by  the  recent  developments  in  the  numerical 
modeling  of  transient  nonlinear  dynamics  phenomena  such  as 
those  accompanying  ballistic  and  blast-loading  conditions. 
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However,  the  tools  used  in  these  analyses  themselves  suffer 
from  a  number  of  deficiencies/limitations  which  prevent  these 
analyses  from  being  more  widely  utilized.  One  of  the  main 
deficiencies  stems  from  the  inability  of  currently  available 
material  models  to  realistically  represent  the  response  of  these 
materials  under  high-deformation  rate,  large-strain,  and  high- 
pressure  loading  conditions,  the  conditions  typically  encoun¬ 
tered  during  projectile  impact  events.  One  of  the  reasons  for  the 
indicated  shortcomings  of  the  (continuum-type)  material  mod¬ 
els  currently  in  use  is  a  lack  of  inclusion  of  the  contributions  of 
various  phenomena  and  processes  occurring  at  different  length 
scales  to  the  overall  behavior  of  the  material  under  dynamic 
loading  conditions.  In  the  present  study,  an  attempt  is  made  to 
upgrade  an  existing  high-rate  composite  material  model 
through  incorporation  of  the  fine-scale  microstructural  (includ¬ 
ing  defect)  effects.  Hence,  the  two  main  topics  covered  in  the 
present  study  include  (a)  continuum-level  composite  material 
models;  and  (b)  multi-length  scale  character  of  the  composite- 
material  microstructure.  These  topics  are  briefly  reviewed  in  the 
remainder  of  this  section. 

1.1  Continuum-Level  Composite  Materia!  Models 

A  review  of  the  open  literature,  carried  out  as  part  of  the 
present  study,  identified  number  of  continuum-level  material 
models  for  the  composite  materials  under  investigation.  Some 
of  these  material  models  (e.g.,  Ref  1)  focus  on  the  role  of 
various  fiber/matrix  damage  processes  solely  due  to  in-plane 
loading.  Since  under  common  ballistic/blast-impact  loading 
conditions,  through-the-thickness  normal  and  transverse  shear 
stresses/strains  are  expected  to  play  an  important  role,  these 
types  of  models  are  of  little  utility  in  the  computational  analysis 
of  protective  structures.  Among  the  material  models  which 


Journal  of  Materials  Engineering  and  Performance 


Volume  22(3)  March  2013—681 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

MAR  2013  2' REPORT  TYPE 

3.  DATES  COVERED 

00-00-2013  to  00-00-2013 

4.  TITLE  AND  SUBTITLE 

Multi-Length  Scale-Enriched  Continuum-Level  Material  Model  for 
Kevlar-Fiber-Reinforced  Polymer-Matrix  Composites 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS (ES) 

Clemson  University , Department  of  Mechanical  Engineering, 241 
Engineering  Innovation  Building, Clemson, SC, 29634 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS  (ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

ARSTRATT 

1 8 .  NUMBER  1 9a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

15 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


include  the  contribution  of  all  the  attendant  stress/strain 
components  on  the  material  response,  most  fall  into  the 
category  of  continuum  damage  mechanics  models.  In  other 
words,  damage  is  treated  as  a  continuum-filed  variable  which 
evolves  during  loading  and  degrades  locally  the  material’s 
stiffness.  Among  these,  the  models  which  appear  most 
noteworthy  are  (a)  a  plane-stress  model  for  the  composite 
materials  reinforced  with  unidirectionally  oriented  fibers 
(Ref  2);  (b)  a  fully  three-dimensional  composite-material  model 
which  accounts  for  the  progressive  character  of  the  main  fiber- 
damage  modes,  i.e.,  fiber  fracture  under  axial  tension,  fiber 
buckling  under  compressive  loading,  fiber  cutting  under 
transverse  shear,  and  fiber  crushing  under  compressive 
through-the-thickness  strains  (Ref  3,  4);  and  (c)  more  recent 
three-dimensional  (3D)  composite  material  models  which 
include  various  continuously  distributed  static  intralamina  and 
interlamina  damage  mechanisms  (e.g.,  fiber  breakage  within  the 
yams,  fiber/matrix  de-bonding,  diffuse  delamination/interlam¬ 
ina  separation,  etc.)  as  well  as  discrete  damage  modes 
(e.g.,  transverse  micro-cracking)  (Ref  5-7). 

The  composite-material  models  briefly  overviewed  above  all 
suffer  from  a  common  deficiency,  i.e.,  they  fail  to  account  for 
the  effects  of  high-loading  rates  under  ballistic/blast  impact 
conditions.  This  deficiency  is  removed  in  the  composite 
material  model  proposed  by  Yen  (Ref  8)  when  the  contribution 
of  the  entire  strain/stress  field  to  all  the  main  fiber  and  matrix 
failure  modes  is  taken  into  account.  This  model  is  adopted  in 
the  present  study,  and  an  attempt  is  made  to  upgrade  this  model 
to  include  the  effects  of  fine  scale  microstructure  to  the  material 
response  under  dynamic  loading  conditions.  Since  this  model  is 
overviewed  in  great  detail  in  section  2,  it  will  not  be  discussed 
any  further  here. 

1.2  Hierarchical/Multi-Length  Scale  Material  Microstructure 

Filament-based  materials/structures  under  investigation  are 
quite  complex,  and  this  complexity  can  be  linked  to  the 
following  main  sources:  (a)  they  contain  a  hierarchical/multi¬ 
length  scale  architecture;  (b)  their  mechanical  response  is  often 
quite  nonlinear  and  rate/time-dependent;  and  (c)  they  are 
associated  with  a  multiplicity  of  complex  phenomena/processes 
(e.g.,  filament  twisting,  interfilament  friction  and  sliding,  etc.). 
In  our  recent  studies  (Ref  9,  10),  an  attempt  was  made  to  help 
us  clarify  the  nature  of  the  multi-length  scale  hierarchy  of  the 
fabric  structures  or  fiber-reinforced  composite  materials  under 
investigation.  A  set  of  eight  length  scales  (described  later)  was 
established  and  briefly  described.  Schematics  and  explanations 
of  the  set  of  eight  microstructural  length  scales  is  provided  in 
Fig.  1.  The  first  column  in  this  figure  shows  a  set  of  simple 
schematics  of  the  material  micro  structure/architecture  at  a  given 
length  scale  along  with  the  labels  used  to  denote  the  main 
microstructural  constituents.  In  the  second  column,  a  brief 
description  is  provided  of  the  material  model  used  to  capture 
the  material  behavior  at  the  length  scale  in  question.  The  key 
features  of  the  material  microstructure/architecture  and  the 
corresponding  material  model  at  each  of  the  length  scales  listed 
in  Fig.  1  can  be  summarized  as  follows:  (a)  Laminate  length- 
scale :  At  this  length  scale,  the  material  possesses  no  discemable 
microstructural  features,  i.e.,  it  is  completely  homogenized.  An 
example  of  the  laminate-length  scale  composite  material  model 
can  be  found  in  Ref  1 1 ;  (b)  Stacked-lamina  length-scale :  At  this 
length  scale,  the  presence  of  discretely  stacked  laminae  is 
recognized,  while  the  material  within  each  lamina  as  well  as 


interlamina  boundaries  are  kept  featureless/homogenized.  An 
example  of  the  composite  material  model  at  this  length  scale 
can  be  found  in  Ref  7;  (c)  Single-lamina  length-scale'.  In  this 
case,  the  micro  structure/architecture  of  the  reinforcement  and 
matrix  phases  are  explicitly  taken  into  account,  while  the  two 
associated  materials  are  treated  as  featureless/homogenized. 
The  resulting  “ two-phase ”  composite  lamina  model  for  each 
lamina  is  then  combined  with  a  fully  homogenized  lamina/ 
lamina  interface  model  to  form  a  stackedlamina  composite 
laminate  structure.  An  example  of  this  type  of  material  model 
can  be  found  in  Ref  12;  (d)  Fabric  unit-cell  length-scale'.  At 
this  length  scale,  a  closer  look  is  given  to  the  architecture  of  the 
woven  fabric.  Specifically,  details  of  yam  weaving  and 
crimping,  yam  cross-section  change,  and  yam  sliding  at  the 
warp-yam/weft-yam  cross-over  points  are  taken  into  account. 
Each  yam  is  fully  homogenized  in  this  case  as  are  the  matrix 
and  the  lamina/lamina  interfaces.  An  example  of  the  model  of 
this  type  can  be  found  in  Ref  13-15;  (e)  Yarn  length-scale :  In 
this  case,  the  internal  stmcture/architecture  of  each  yam  is 
accounted  for  explicitly.  In  other  words,  yams  are  considered  as 
assemblies  of  nearly  parallel  fibers/filaments  which  are 
mechanically  engaged  by  either  the  application  of  a  light  twist 
to  the  yam  or  by  wrapping  a  fiber  around  the  fiber/filament 
assembly.  Typically,  the  detailed  microstmcture  of  the  yams  is 
not  incorporated  into  material  models  used  in  large-scale 
simulations  of  projectile/protection- system  impact  interactions 
(due  to  the  unmanageably  large  computer  resources  needed). 
Instead,  computational  results  of  the  mechanical  response  of  the 
individual  yams,  when  subjected  to  a  variety  of  loading 
conditions,  are  used  to  enrich  the  homogenized  material  model 
for  the  yams.  The  latter  is  subsequently  used  at  the  fabric  unit¬ 
cell  length  scale.  An  example  of  this  approach  can  be  found  in 
Ref  16;  (f)  Fiber  length-scale'.  At  this  length  scale,  fibers  are 
considered  as  assemblies  of  aligned  long-chain  molecules 
which  are  held  together  by  non-bonding  (van  der  Waals  or 
Coulomb)  forces.  In  the  case  of  Kevlar®  type  fibers,  there  is  a 
substantial  experimental  support  (e.g.,  Ref  17)  for  the  existence 
of  fibrils  within  the  fibers.  Fibrils  are  smaller  bundles  of 
molecular  chains  within  which  chain  molecules  are  tightly 
bonded  into  a  perfect  or  nearly  perfect  crystalline  phase.  Thus, 
in  the  case  of  Kevlar®  type  fibers,  fibers  can  be  considered  as 
an  assembly  of  fibrils;  (g)  Fibril-level  length-scale'.  As 
mentioned  above,  within  the  fibrils,  the  molecular  structure  is 
crystalline  or  nearly  crystalline.  However,  as  clearly  demon¬ 
strated  in  our  recent  study  (Ref  9),  the  material  at  this  length 
scale  (as  well  as  the  fiber  length  scale)  may  contain  a  variety  of 
micro  structural  and  topological  defects  and  chemical  impurities 
which  may  significantly  alter  its  properties.  As  in  the  fiber 
length  scale  case,  the  material  is  treated  as  a  collection  of 
discrete  interacting/bonded  particles  and  analyzed  using  the 
atomic/molecular  modeling  tools/procedures;  and  (h)  Molecu¬ 
lar  chain-level  length-scale'.  At  this  length  scale,  chemical 
stmcture  and  conformation  of  the  individual  molecules  consti¬ 
tuting  the  chain  are  analyzed  using  the  aforementioned  atomic/ 
molecular  modeling  tools/procedures.  The  main  goal  of  the 
analysis  at  this  length  scale  is  to  identify  the  most-likely 
molecular  conformations  present  in  the  fibrils. 

1.3  Objective 

As  mentioned  earlier,  the  main  objective  of  the  present  study 
is  to  identify  and  quantify  fine-scale  microstructural  effects  on 
the  dynamic  mechanical  response  of  the  composite  materials 
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Fig.  1  Various  length  scales  and  the  associated  material  model  assumptions/simplifications  used  in  the  study  of  polymer-matrix  composite 
materials  with  high-performance  fiber-based  structures 


under  investigation  and  to  incorporate  these  effects  into  the 
continuum-level  material  model  of  Yen  (Ref  8). 

1.4  Organization  of  the  Paper 

A  brief  overview  of  the  continuum-level  composite  material 
model  proposed  by  Yen  (Ref  8)  is  provided  in  section  2.  A 
summary  of  the  fiber  and  sub-fiber  level  microstructure  and 
morphology  of  Kevlar® -based  materials  including  the  polymer- 
synthesis  and  fiber  fabrication  processes-induced  defects  is 
presented  in  section  3.  The  procedure  used  to  upgrade  the 
model  of  Yen  (Ref  8)  and  the  results  obtained  using  the 
upgraded  material  model  are  all  presented  and  discussed  in 
section  4.  The  key  findings  resulting  from  the  present  study  are 
summarized  in  section  5. 


2.  Overview  of  the  Yen  Composite  Material 
Model  (Ref  8) 

In  this  section,  a  brief  overview  is  provided  of  the  rate- 
dependent  composite  material  model  developed  by  Yen  (Ref  8). 
The  defining  features  of  this  model  can  be  summarized  as 
follows: 

(a)  The  model  is  an  extension  of  the  Hashin  composite- 
material  failure  model  (Ref  18)  which  more  accurately 
accounts  for  the  three-dimensional  nature  of  the  stress/ 
strain  states,  high  strain  rates,  and  high  pressures 
encountered  during  ballistic  impact  of  composite  lami¬ 
nates  reinforced  with  plain-woven  fabric; 
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Yarn  Length-scale 

Each  yam  is  considered  as  an  assembly 
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engaged  discrete  fibers 
Fiber-fiber  contact  and  sliding  are 
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Fiber  Length-scale 

Fibers  are  considered  as  an  assembly  of 
n early-parallel  fibrils  held  together  by 
weak  van  der  Walls  forces 
Adjacent  fibrils  differ  with  respect  to  the 
orientation  of  the  stacked  sheets. 


Fibril  Length-scale 

Fibrils  are  considered  as  crystalline 
thread-like  entities  containing  parallel 
sheets  of  aligned  and  hydrogen 
bonded  molecular  chains. 

Inter-sheet  bonding  is  of  a  van  der  Walls 
character. 


Molecular-level  Length-scale 

•  Individual  molecules  are  represented  as 

assemblies  of  covalently  bonded 
atoms. 

*  The  effect  of  hydrogen  bonding  on  the 

molecular  topology  (e  g.  the  formation 
of  helical  structure  is  considered 
explicitly). 


Fig.  1  continued 


(b)  The  material  model  is  strain  based  and  utilizes  the  con¬ 
cept  of  a  ply/lamina-averaged  local  strain; 

(c)  The  model  incorporates  Mohr-Coulomb  (Ref  19)  shear 
strength  effects,  i.e.,  it  accounts  for  the  fact  that  the 
(matrix)  in-plane  shear  strength  is  generally  a  function 
of  the  (tensile/compressive)  character  and  magnitude  of 
the  through-the-thickness  lamina  strain; 

(d)  Within  the  model,  no  distinction  is  made  between  fibers 
and  yams.  Hence,  to  be  consistent  with  the  original 
model,  in  this  section,  the  term  “fiber”  will  be  used  to 
denote  the  thread-like  entity  which  is  used  in  the  fabric¬ 
weaving  process.  In  the  next  and  the  subsequent  sec¬ 
tions,  the  same  entity  will  be  referred  to  as  a  “yarn”; 

(e)  The  following  three  fiber  and  two  matrix  failure  modes 
are  considered:  (i)  fiber  failure  due  to  a  combination  of 
axial  tensile  and  transverse  shear  strains;  (ii)  fiber  buck¬ 
ling  failure  due  to  axial  compressive  strains;  (iii)  fiber 
cmshing  failure  due  to  lamina  through-the-thickness 
compressive  strains;  (iv)  matrix  in-plane  shear  failure; 


and  (v)  matrix  delamination  failure  due  to  a  combined 
effect  of  the  through-the-thickness  tensile  strains  and 
transverse  shear  strains;  and 

(f)  The  model  is  essentially  of  a  linear-elastic  orthotropic 
type  with  progressively  degradable  elastic  stiffness  con¬ 
stants.  In  other  words,  a  generalized  Hooke’s  law  relation 
is  assumed  to  hold  between  stress  and  strain  quantities, 
while  damage  is  treated  as  a  process  which  takes  place 
gradually  during  ballistic  impact.  Consequently,  the  mod¬ 
el  relies  on  the  definition  of  the  damage  initiation  criteria 
and  the  corresponding  damage  evolution  laws. 

2.1  Damage  Initiation  Criteria 

In  this  section,  a  brief  overview  is  provided  of  the  functional 
relations  defining  the  onset  (and  continuation)  of  damage  for 
each  of  the  five  aforementioned  fiber/matrix  failure  modes.  For 
each  of  the  failure  modes  in  question,  as  will  be  shown  shortly, 
because  distinction  is  made  between  weft  and  warp  fibers,  a  set 
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of  seven  equations  had  to  be  used  to  fully  describe  damage 
initiation.  Each  of  the  seven  equations  has  a  general  form: 
f  —  rf  =  0  (7=1,..., 7)  where  /  represents  the  appropriate 
strain-based  loading  function,  while  rt  represents  the  corre¬ 
sponding  strain-based  damage  threshold.  Initially,  each  rt  is  set 
to  a  value  of  1.0  and,  as  a  result  of  progressive  damage,  rt 
values  continue  to  increase. 

It  should  be  noted  that  in  the  remainder  of  this  article,  the 
following  Cartesian  coordinate  system  is  used:  x  and  y  are, 
respectively,  aligned  with  the  weft  and  the  warp  fibers,  while  z  is 
aligned  in  the  composite- laminate  through-the-thickness  direction. 

2.1.1  Fiber  Damage  Under  Combined  Axial  Tension  and 
Transverse  Shear.  The  onset  and  progression  of  this  failure 
mode  is  assumed  to  be  defined  by  a  superposition  of  the 
quadratic  contributions  of  the  axial  tensile  and  transverse  shear 
strains.  For  the  weft  and  the  warp  fiber  directions,  the 
corresponding  damage  initiation  functions  can  be  written  as 

/,->-?=  (^)+(^)2-^  =0  (Eq  1) 


fi-rl 


^y{Zy)\  .  fGyz£yz\  „2  —  0 


(Eq  2) 


where  E  denotes  the  fiber  Young’s  modulus,  G  the  corre¬ 
sponding  shear  modulus,  S  the  fiber  strength,  8  the  relevant 
strain  component,  the  Macaulay  brackets  ()  are  used  to  de¬ 
note  the  ply/lamina-average  local  quantity,  and  subscripts  T 
and  FS  denote,  respectively,  tension  and  fiber  shear.  It  should 
be  noted  that  in  Eq  1  and  2,  the  laminate  transverse  shear 
strength  is  assumed  to  be  controlled  by  the  corresponding 
shear  strength  of  the  fiber  itself. 

2.1.2  Fiber  Damage  Under  Axial  Compression.  The 
onset  and  progression  of  this  failure  mode  is  assumed  to  be 
defined  by  the  compressive  axial  strains,  effect  of  which  is 
modulated  by  the  lamina  through-the-thickness  compressive 
strain.  In  other  words,  axial  compression-induced  fiber  buck¬ 
ling  is  resisted  by  the  stabilizing  effect  of  the  through-thickness 
compressive  strain.  For  the  weft  and  the  warp  fiber  directions, 
the  corresponding  damage  initiation  functions  can  be  written  as 


I2 

Sxc 


—r\  —  0 


(Eq  3) 


f\  -r\  = 


(Eq  4) 


where  prime  is  used  to  denote  modulation  of  the  axial  com¬ 
pressive  strains  as  =  —  ex  —  (—8 Z){EZ/EX)  and  s'y  =  — Sy— 
(—8 Z)(EZ/Ey),  and  the  subscript  C  denotes  compression. 

2.1.3  Fiber  Damage  Under  Transverse  Compressive 
Loading.  This  failure  mode  involves  fiber  crushing  under 
sufficiently  high  through-the-thickness  compressive  strain,  and 
its  onset/progression  is  governed  by  the  following  equation: 

/5-^=(F^)2-r^  =  0  (Eq  5) 

where  the  subscript  FC  denotes  a  fiber  compression  quantity. 

2.1.4  Matrix  In-plane  Shear  Failure.  This  mode  of 
failure  is  associated  with  the  fact  that  the  matrix  shear  strength 
is  lower  than  the  fiber  shear  strength  and  that  the  composite 


material  in  question  possesses  a  laminated  architecture.  Con¬ 
sequently,  sufficiently  high  in-plane  shear  strains  may  cause 
matrix  shear  failure  without  giving  rise  to  any  significant  fiber 
damage.  The  onset  and  progression  of  this  type  of  damage  are 
governed  by  the  following  functional  relationship: 

/6-^(^)2-^0  (Eq  6) 


2.1.5  Matrix  Delamination  due  to  Through-the-Thick- 
ness  Tensile  and  Shear  Stresses.  This  mode  of  failure  also 
gives  rise  to  the  matrix  in-plane  failure,  and  its  onset/ 
progression  is  described  by  the  following  relation: 


fi-t=S2 


%{ez)\ 2  (  GyzSyz  V  /  GzxZzx  \ 2 1 
,  SzT  )  \SyzO  +  ^SR  /  \Szx  0  +  ^SR  /  J 


-ri  =  0 


(Eq  7) 


Examination  of  Eq  7  reveals  the  presence  of  three  terms  within 
the  curly  brackets  on  its  right-hand  side.  The  first  term  defines  the 
contribution  of  tensile  through-the-thickness  strains,  while  the 
remaining  terms  represent  the  contributions  of  the  transverse  shear 
strains,  to  matrix  delamination.  In  Eq  7,  subscript  0  is  used  to 
denote  the  appropriate  shear  strength  in  the  absence  of  the 
associated  compressive  normal  strain,  while  subscript  SR  is  used  to 
denote  an  increase  in  this  shear  strength  due  to  the  presence  of  a 
normal  compressive  strain.  The  SSr  (sliding  resistance  term  is 
assessed  using  the  Mohr-Coulomb  approach  as 

Ssr  =  Ez  tan  4>(-sz)  (Eq  8) 

where  tanc|)  is  the  friction  angle. 

Equation  7  also  contains  a  stress-concentration  (scaling) 
factor,  S.  This  factor  is  assigned  a  value  larger  than  1.0  in  the 
non-delaminated  portion  of  the  material  adjacent  to  the 
delamination  boundary/front.  This  algorithm  has  been  found 
to  introduce  enough  flexibility  into  the  model  so  that  a  better 
agreement  can  be  obtained  with  the  corresponding  experimen¬ 
tal  results  pertaining  to  the  temporal  evolution  of  the  delam¬ 
ination  zone  (Ref  8); 


2.2  Damage  Evolution  and  Material  Degradation 

2.2.1  Damage  Surface.  Equations  1-7  given  above 
define  within  the  strain  space  a  damage  initiation/progression 
surface.  Examination  of  these  equations  reveals  that  the 
material  stiffness  is  characterized  by  three  Young’s  moduli 
and  three  shear  moduli.  The  initial  damage  initiation/progres¬ 
sion  surface  is  defined  by  substituting  the  initial  values  of  the 
six  moduli  into  Eq  1-7.  As  the  material  moduli  degrade,  the 
failure  surface  within  the  strain  space  expands  (i.e.,  additional 
straining  is  required  to  produce  additional  damage).  Within  the 
composite  material  model  developed  by  Yen  (Ref  8),  damage 
progression  is  assumed  to  cause  degradation  in  material 
stiffness  (without  an  accompanying  loss  in  material  strength) 
properties.  The  same  assumption  will  be  adopted  in  the  present 
study.  In  our  forthcoming  study,  however,  it  is  planned  to 
address  the  problem  of  strength  degradation.  It  should  be  noted 
that  although  stiffness  is  considered  as  a  non-degradable 
material  property,  it  will  be  considered  as  a  stochastic  rather 
than  a  deterministic  quantity  in  the  present  study. 

2.2.2  Stiffness  Degradation  Variables.  To  account  for 
the  aforementioned  damage-induced  material- stiffness  degra¬ 
dation,  six  stiffness  degradation  variables,  coj  (J  =  1,...,6),  are 
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first  introduced.  Then,  the  effect  of  damage  on  the  material 
elastic  compliance  matrix  is  described  by  the  following  relation 
(Ref  2): 
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(Eq9) 


Finally,  the  “ degraded ”  material  elastic  stiffness  matrix  is 
obtained  by  inverting  the  material  elastic  compliance  matrix  as 
defined  by  Eq  9.  This  procedure  shows  that  to  compute  the 
instantaneous  (degraded)  material  stiffness  properties,  one  must 
know  the  corresponding  values  of  the  six  stiffness  degradation 
variables.  To  evaluate  the  same,  one  must  define  the  appropriate 
evolution  equations  for  the  stiffness  degradation  variables. 

Following  Matzenmiller,  et  al.  (Ref  2),  the  needed  evolution 
equations  are  defined  as 

Y  °'9/<  (Eq  !0) 

i 

where  (j),  (i=  1,...,7)  denotes  the  rate  of  progression  of  each 
of  the  failure  modes  (the  initiation  of  which  was  described  by 
Eq  1-7),  while  qjt  represents  a  binary  (0  or  1)  two  dimen¬ 
sional  array  with  j  =  1,...,6  and  quantifies  the  contribution  of 
different  failure  modes  to  the  evolution  of  stiffness  degrada¬ 
tion  variables. 

2.2.3  Damage  Evolution  Laws.  Next,  cj)z  functional  rela¬ 
tions  are  derived  for  each  of  the  failure  modes.  It  should  be  first 
noted  that  failure  and  the  associated  stiffness  degradation  are 
irreversible  and  that,  during  loading,  an  increase  in  the  extent  of 
damage  of  a  given  mode  occurs  only  if  the  resulting  strain 
increment  crosses  the  corresponding  failure  surface,  i.e.,  that 
the  projection  of  the  strain  increment  onto  the  local  failure 
surface  outward  normal  is  positive  as  0, 

where  k  (=  1,...,6)  is  used  to  denote  six  components  of  the 
strain  vector.  Then,  by  defining  the  failure-mode-specific 
progression  function,  yf,  the  rate  of  progression  of  failure 
mode  i  can  be  defined  as 


(Eq  11) 


where  no  summation  is  implied  over  repeating  indices.  By 
choosing 

Y,=^(l-rf_1  (Eq  12) 

where  the  material  parameter  m  quantifies  the  stiffness  sensi¬ 
tivity  to  damage  extent  and  noting  that  ^ k  (dfi/dzk)Zk  =  f 
the  final  expression  for  the  rate  of  failure  progression  can  be 
defined  as 


(j)z-  =  -  (1  —  lfi  (no  summation  over /)  (Eq  13) 

2.2.4  Derivation  of  the  Damage  Coupling  Matrix.  To 

assemble  the  damage  coupling  matrix,  (j  denotes  the  elastic 
moduli,  while  i  denotes  the  failure  mode),  the  following 
approach  was  taken: 


(a)  Failure  of  the  weft  and  warp  fibers  both  in  tension  and 
transverse  shear  and  in  compression  degrades  the  associ¬ 
ated  Young’s  modulus  and  shear  moduli; 

(b)  Fiber  crushing  mode  degrades  all  six  elastic  moduli; 

(c)  In  the  case  of  the  in-plane  shear  failure,  only  the  associ¬ 
ated  shear  modulus  is  affected;  and 

(d)  Delamination  degrades  the  through- the-thickness  Young’s 
modulus  as  well  as  the  two  transverse  shear  moduli. 


The  resulting  damage  coupling  matrix  is 


0  1  0 
1  0  1 
0  0  0 
1  1  1 
1  0  1 
0  1  0 


1  0  0 
1  0  0 
1  0  1 
1  1  0 
1  0  1 
1  0  1 


(Eq  14) 


2.2.5  Integrated  Form  of  the  Stiffness  Degradation 
Variables.  By  combining  Eq  10,  12,  and  13,  the  following 
expression  is  obtained  for  the  contribution  of  the  failure  mode, 
i,  to  each  stiffness  degradation  variable: 

coy  =  1 -e»(‘ "'I"),  n>\  (Eq  15) 

Examination  of  Eq  15  reveals  that  initially,  when  rt  is  equal 
to  1.0,  the  corresponding  contribution  to  coj  is  zero.  During 
loading,  when  the  strain  increment  causes  the  strain  to  fall 
outside  the  damage  surface  associated  with  failure  mode  /,  the 
corresponding  damage  surface  crossed  is  expanded  outward. 
This,  in  accordance  with  Eq  1-7,  gives  rise  to  an  increase  in 
and,  in  turn,  via  Eq  15  to  an  increase  in  coj. 

2.2.6  Residual  Stiffness.  According  to  Eq  15,  at  suffi¬ 
ciently  high  rb  the  stiffness  degradation  variable  may  acquire  a 
value  of  1.0  in  which  case  the  material  completely  loses  the 
associated  load  carrying  capacity.  While  this  is  physically 
realistic  for  some  of  the  failure  modes,  e.g.,  the  fiber  failure 
under  axial  tension/transverse  shear,  it  may  not  be  fully  realistic 
under  other  failure  modes,  e.g.,  fiber  axial  compression. 
Consequently,  the  following  modifications  are  made  in  Eq  15: 


(a)  For  the  fiber  axial  compression  failure  mode: 
CD/  =  ^(1_rr)^ (i  _  r|.),  z  =  3or4 


(Eq  16) 


where  the  subscript  R  is  used  to  denote  the  respective  resid¬ 
ual-strength  quantities;  and 


(b)  In  the  case  of  delamination  failure  mode,  distinction  is 
made  between  this  failure  mode  under  compressive  and 
tensile  through- the-thickness  strains.  In  the  former  case, 
delamination  surfaces  are  assumed  to  be  closed,  and 
Eq  15  for  the  two  transverse  shear  moduli  is  modified 
to  include  the  effect  of  the  sliding  resistance  as 

_ ^sr _ 

[{Gyzeyz)2+(Gzxezx)2Y2  j 

(Eq  17) 
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.  On  the  other  hand,  in  the  case  of  the  delamination  under 
tensile  through-the-thickness  strains,  no  modifications  are  made 
to  Eq  15. 


2.3  Rate  Dependency  of  Strength  and  Stiffness 


As  mentioned  earlier,  one  of  the  defining  features  of  the 
model  developed  by  Yen  (Ref  8)  is  the  inclusion  of  material 
(strength  and  stiffness)  rate  dependency.  To  account  for  rate 
dependency  of  different  components  of  material  strength,  the 
following  scaling  relation  was  used: 

{SRT}  =  {So}  (l  +  Cl  InEf)  (Eq  18) 

where 
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Ci  is  the  strength  rate  constant,  and  the  6  x  1  column  vector 
{S0}  is  equal  to  {SKT}  when  each  of  the  relevant  strain  rates 
is  equal  to  the  single  reference  strain  rate  So  =  1.0s-1. 

The  rate  dependency  of  different  stiffness  components  was 
accounted  for  in  a  similar  fashion  as 

{£„}  =  {£0}  (l  +  C2  lnEf)  (Eq  19) 

where 
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C2  is  the  stiffness  rate  constant,  and  the  6  x  1  column  vector 
{E’o}  is  equal  to  {Ert}  when  each  of  the  relevant  strain  rates 
is  equal  to  the  single  reference  strain  rate  8q  =  1.0s-1. 


2.4  Material  Model  Parameterization 

The  material  model  presented  above  contains  14  strength 
parameters  (SxT,  SyT,  SxF s,  SyF5,  *S*o  Syo  SFC,  Sxy  SzF,  Syzo>  Szxo, 
phi,  SxCr,  and  SycR),  six  stiffness  parameters  (Ex,  Ey,  Ez,  Gxy, 
Gyz ,  and  Gzx);  two  rate  dependency  parameters  (Q  and  C2),  and 
two  damage  parameters  (S  and  m).  These  parameters  are 
commonly  assessed  by  applying  regression  analysis  to  the 
results  from  a  variety  of  mechanical  tests  over  a  range  of  strain 
rates.  The  large  number  of  the  model  parameters  makes  this 
parameterization  effort  quite  cumbersome,  time  consuming  and 
costly.  In  addition,  as  will  be  discussed  in  the  next  section, 
some  of  these  parameters  are  stochastic  quantities  and  their 
complete  specification  requires  the  use  of  the  appropriate 
probability  distribution  functions.  Determination  of  these 
distribution  functions  is  even  more  costly.  Hence,  a  multi¬ 
length  scale  computational  approach  is  introduced  in  the 
next  section  which  allows  establishment  of  the  appropriate 


distribution  functions  through  the  use  of  time-efficient,  low  cost 
computational  methods  and  tools. 


3.  Microstructure/Property  Relations  in  Kevlar® 
Fibers 

In  this  section,  a  brief  overview  is  provided  of  the  basic 
micro  structural  features  in  Kevlar®  at  the  fibril/fiber/yam  length 
scale(s),  and  their  effects  on  the  key  material  properties. 

3.1  Fiber  Crystalline  Structure 

Kevlar®  (p-phenylene  terephthalamide,  PPTA)  belongs  to 
the  family  of  polymeric  materials  known  as  polyamides. 
Polyamides  are  typically  classified  as  aromatic  polyamides  or 
aramids  (e.g.,  Kevlar®,  the  subject  of  the  current  study)  and 
non-aromatic  polyamides  (e.g.,  Nylon  6,6).  Kevlar®  chains  are 
inherently  stiff  and  hence  the  chains  do  not  readily  flex. 
Consequently,  in  contrast  to  the  flexible  polymeric  molecules 
which  can  undergo  extensive  folding  and  form  the  common 
(crystalline  +  amorphous)  two  phase  structure,  the  Kevlar® 
fibers  tend  to  typically  form  either  a  paracrystalline  or  a  fully 
crystalline  structure.  The  tendency  for  the  formation  of 
paracrystalline  or  crystalline  structures  is  promoted  by  the 
presence  of  the  planar  phenylene  and  amide  groups  and  by 
the  ability  of  the  adjacent  chains  to  form  hydrogen  bonds.  In  the 
case  of  the  paracrystalline  structure,  Kevlar®  molecules  are  all 
aligned  in  the  same  direction  but  no  order  exists  in  a  plane 
orthogonal  to  this  direction.  In  sharp  contrast,  in  the  case  of  the 
fully  crystalline  Kevlar®  fibers,  molecules  are  aligned  in  all 
three  mutually  orthogonal  directions. 

The  Kevlar®  crystal  structure  is  layered  and  consists  of 
parallel  (ABABAB...)  stacked  sheets  formed  due  to  the 
hydrogen  bonding  between  the  adjacent  parallel  Kevlar® 
molecules.  The  intersheet  bonding  on  the  other  hand  is  mainly 
of  the  van  der  Waals  (and  p-electron  weak  chemical-bond)  type. 
Owing  to  the  relatively  low  strength  of  the  van  der  Waals  and 
the  p-electron  type  weak  intersheet  bonding,  Kevlar®  fibers  are 
prone  to  the  formation  of  stacking  faults  and  kink  bands,  and 
typically  possess  an  inferior  compressive  strength  and  buckling 
resistance. 

3.2  Common  Defects  in  Kevlar*  Fibers 

As  in  most  engineering  materials,  properties  of  Kevlar® 
fibers  are  greatly  affected  by  the  presence  of  various  defects/ 
flaws.  These  flaws  can  be  introduced  during  Kevlar®  synthesis 
and  Kevlar®  fiber  fabrication  processes.  In  our  recent  study 
(Ref  10),  a  detailed  overview  was  provided  of  the  main  steps 
associated  with  Kevlar®  synthesis  and  Kevlar®  fiber  process¬ 
ing.  In  the  next  section,  a  brief  summary  is  provided  of  the 
commonly  observed  synthesis/processing-induced  defects  in 
Kevlar®  fibers. 

A  summary  of  the  Kevlar®  fiber’s  most  common  defects, 
their  dimensionality,  their  cause,  methods  of  reducing  their 
density,  and  their  typical  concentrations  is  provided  in  Table  1 . 
Molecular-level  schematics  of  these  defects  are  displayed  in 
Fig.  2(a)-(h).  These  defects  are  expected  to  have  a  profound 
effect  on  the  fiber  properties,  as  well  as  on  the  properties  of 
coarser-scale  materials  and  structures  containing  Kevlar® 
fibers  (e.g.,  yams,  fabrics,  plies,  lamina,  and  laminates). 
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Table  1  A  summary  of  the  most  common  defects  found  in  Kevlar®-based  materials 


Flaw  group 

Flaw 

Causes 

How  causes  can  be  remedied/ 
mitigated 

Range  for  concentration 

Isolated  chain 

-COOH 

H2S04  catalyzed  hydrolysis 

Use  concentrated  H2S04  for 

0.35  per  PPTA  chain 

ends  (point 

causing  PPTA  chain  scission. 

dope  preparation.  Shorten  the 

for  each  flaw(a) 

defect) 

Na+  deficiency  with  respect 
to  complete  neutralization  of 
side/end  acidic  groups 

fiber  wash  time 

(~  350  ppm-mass- 
based) 

-nh2 

H2S04  catalyzed  hydrolysis 
causing  PPTA  chain  scission. 
Na+  deficiency  with  respect 
to  complete  neutralization  of 
side/end  acidic  groups 

Use  higher  concentration 

NaOH  solution 

0.35  per  PPTA  chain 
for  each  flaw(a) 

(~  350  ppm-mass- 
based) 

-COO“Na+ 

COOH  neutralization  with  Na+ 

No  remedy  required  since 
this  is  one  of  the  preferred 
chain  ends 

1 . 1  per  PPTA  chain(a) 

(~  1100  ppm-mass- 
based) 

-NH3+HS04“ 

Sulfonation  of  the  NH2  chain 
ends 

Increase  the  H2S04  removal 
and  neutralization  rate 

0.2  per  PPTA  chain(a) 

(~  200  ppm-mass- 
based) 

Side  groups 

-so3h 

Exposure  of  PPTA  in  the  dope 

Reduce  the  H2S04  concen¬ 

~1300  ppm 

(point 

defect) 

to  concentrated  H2S04 
(sulfonation) 

tration  in  the  dope 

(mass-based) 

-S03“Na+ 

Neutralization  of  sulfonic  acid 
side  groups  by  NaOH 

Remedy  may  not  be  required 
since  this  side  group 
improves  fiber  longevity. 
However  mechanical  perfor¬ 
mance  may  be  compromised 

~2500  ppm 
(mass-based) 

Voids  and 

Microvoids 

Swelling  induced  by  hydration 

Increase  the  extent  of  sodium 

~150  ppm 

interstitials 
(point  defects) 

of  intra-fibrillar  Na2S04 

salt  dissolution  by  prolonged 
exposure  of  fibers  to  boiling 
water 

(mass-based) 

Mobile  trapped  H2S04 

Non-neutralized  or  unwashed 
intra-fibrillar  H2S04 

Thorough  washing  in  hot 
solvent  aqueous  bath 

~70  ppm  (mass-based) 

Defect  bands 

NH3+HS04“ 

Coulombic  forces-induced 

The  phenomenon  is  not 

One  band  every 

(planar  defects) 

agglomerated 

clustering  of  ion-terminated 

well  understood  so  no 

40-60  nm  of  fibril 

(a)  Extruded  fibers 

chain  ends 

chain  ends 

remedy  is  obvious 

(ca.  3000  ppm- 
mass-based) 

Consequently,  a  high-fidelity  material  model  for  Kevlar®  fiber- 
reinforced  composite  materials  must  account  for  the  presence, 
distribution,  and  potency  of  these  defects  in  modifying  material 
response  under  dynamic  loading  conditions.  As  mentioned 
earlier,  the  main  objective  of  the  present  study  is  to  upgrade  one 
of  the  currently  used  continuum-level  composite  material 
models  via  incorporation  of  the  effect  of  these  defects. 

3.3  Effect  of  Defects  on  Kevlar"  Properties 

In  this  section,  a  brief  summary  is  provided  of  the  main 
results  obtained  in  our  recent  study  (Ref  10)  in  which 
molecular-level  modeling  was  used  to  establish  the  effects  of 
various  defects  on  the  material  properties  (primarily  strength 
and  stiffness).  In  Ref  10,  a  procedure  was  introduced  to 
establish  relationships  between  topology/kinematics  and  prop¬ 
erties  of  Kevlar®  fibrils,  fibers,  and  yams,  Fig.  3(a)  and  (b).  In 
Fig.  3(a),  a  Kevlar®  fiber  is  depicted,  and  it  is  shown  that  it 
consists  of  a  number  of  nearly  parallel  fibrils.  Within  each  of 
the  fibrils,  the  fibril’s  axis  lies  within  the  Kevlar®  sheets,  and 
the  fibrils  merely  differ  in  the  orientation  of  the  sheet  normals 
(the  sheet  normal  in  each  fibril  is  orthogonal  to  the  fibril  axis). 
Figure  3(b)  displays  the  assumed  relationship  between  the 
fibers  and  the  yams.  Simply  stated,  the  yams  are  treated  as 


assemblies  of  nearly  parallel  fibers  which  are  given  a  light  twist 
to  mechanically  engage  the  fibers. 

Within  the  analysis  presented  in  Ref  10,  it  was  assumed  that 
the  longitudinal  stiffness,  strength,  and  ductility  of  the  fibrils 
make  the  dominant  contribution  to  the  corresponding  fiber’s 
and,  in  turn,  yam’s  properties.  As  far  as  the  transverse  strength 
of  fibers  and  yams  are  concerned,  they  are  expected  to  be 
controlled  by  weak  interfibril  (and  in  the  case  of  yams,  weak 
interfiber)  van  der  Waals/mechanical  interlocking  forces  and 
to  be  less  dependent  on  the  fibril’s  transverse  mechanical 
properties. 

The  results  obtained  in  Ref  10  can  be  summarized  as 
follows:  (a)  fibril’s  stiffness  was  found  to  be  affected  by  both 
the  type  and  concentration  of  the  defects.  An  example  of  the 
typical  results  pertaining  to  the  effect  of  defect  type/concen¬ 
tration  on  the  Kevlar  fibril’s  stiffness  is  depicted  in  Fig.  4(a) 
and  (b).  In  Fig.  4(a),  the  case  of  -COOH  chain  end  defect  is 
considered,  while,  in  Fig.  4(b),  the  case  of  -S03H  side  group 
defect  is  analyzed;  (b)  As  far  as  the  fibril’s  axial  strength  is 
concerned,  it  was  found  to  be  affected  by  defect  type  and  the 
size  and  composition  of  the  largest  defect  cluster.  Examples  of 
typical  results  obtained  are  shown  in  Tables  2  and  3.  The 
clusters  composed  of  defects  of  the  same  type  are  depicted  in 
Table  2,  while  the  effect  of  mixed-type  defect  (chain  end  and 
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Fig.  2  Examples  of  the  molecular-level  defects  in  Kevlar®  fibers:  (a)  a  COOH  chain-end  defect;  (b)  a  NH2  chain-end  defect;  (c)  a  COOH 
chain-end  defect  cluster;  (d)  a  S03H  side-group  defect;  (e)  a  S03H  side-group  defect  cluster;  (f)  a  microvoid  defect;  (g)  a  H2S04  interstitial 
defect;  and  (h)  a  collection  of  chain  ends  forming  a  defect  band 


side  group)  clusters  is  depicted  in  Table  3.  An  examination  of 
the  results  reported  in  Tables  2  and  3  reveals  that  the  extent  of 
the  loss  of  fibril’s  axial  strength  increases  with  the  cluster  size/ 
strength;  (c)  The  aforementioned  finding  that  the  fibril’s  axial 
strength  is  a  function  of  the  largest  defect  cluster  and  its 
composition  indicated  that  fibril’s  strength  is  a  stochastic 
property.  Hence,  in  Ref  10,  results  like  those  displayed  in 
Tables  2  and  3  are  combined  with  the  prototypical  defect 
concentration  densities,  Table  1,  to  derive  the  associated 
Kevlar®  fibril-strength  distribution  function,  Fig.  5(a)  and  (b); 
(d)  In  accordance  with  Fig.  3(a),  Kevlar®  fibers  are  assumed  to 


be  transversely  isotropic  (with  the  fiber  axis  being  the  unique 
material  direction)  although  the  constituent  fibrils  possess 
orthotropic  symmetry.  Consequently,  (i)  fiber  axial  stiffness  and 
strength,  En,  and  ox  are  set  equal  to  their  fibril’s  counterparts; 
(ii)  the  fiber  transverse  normal  stiffness,  E22  =  E33,  is  set  equal 
to  the  more  compliant  ^33  fibril’s  modulus;  (iii)  all  the  three 
shear  moduli  Gi2,  Gi3,  and  G2 3  are  set  equal  to  the  same  value 
which  is  controlled  by  the  compliant  Kevlar®  intersheet  sliding 
resistance  (i.e.,  G12  of  the  fibrils);  and  (e)  To  assess  the 
yam-level  stiffness  and  strength  properties  from  the  known 
fiber  properties,  additional  effects  arising  from  the 
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Fig.  2  continued 

fiber-twist-induced  interfiber  friction  had  to  be  taken  into 
account.  This  procedure  mainly  affected  the  yam  axial  strength 
(or  more  precisely,  yam  axial  strength  distribution).  It  was 
generally  found  that  the  fibril-level  axial  strength  distribution 
function  displayed  in  Fig.  5(a)  and  (b)  was  modified  in  the  yam 
case  by  an  interfiber  friction  dependent  (positive)  strength  shift. 
For  example,  for  a  prototypical  extent  of  fiber  twist,  an  increase 
in  the  strength  of  ca.  0.2  GPa  was  found. 

4.  Enrichment  of  Yen  (Ref  8)  Composite 
Material  Model 

In  this  section,  some  of  the  findings/results  reported  in 
section  3  are  used  to  enrich  the  continuum-level  material  model 
developed  by  Yen  (Ref  8).  Specifically,  Kevlar®-fiber’s  strength 
distribution  functions  like  the  one  displayed  in  Fig.  5(a)  and  (b) 
are  used  to  impart  a  stochastic  character  to  the  Kevlar®-fiber’s 
tensile  strength  and  transverse- shear  strength  parameters  SxT, 
SyT,  SxF s,  and  SyFS.  It  should  be  recalled  that  these  strength 
parameters  are  treated,  both  in  the  material  model  developed  by 
Yen  (Ref  8)  and  in  its  present  modification,  as  non-degradable 
material  properties.  To  demonstrate  the  effect  of  material  model 


Defect  band 


modifications  on  the  material  mechanical  response,  the  case  of 
uniaxial  tension  in  the  weft  direction  is  first  considered.  This  is 
followed  by  a  dynamic  computational  analysis  in  which  a  rigid 
projectile  of  a  right  circular  cylindrical  geometry  is  impacting, 
at  a  zero  obliquity  angle,  a  target  laminate  of  rectangular- 
parallelepiped  geometry.  In  the  latter  case,  it  was  critical  to 
demonstrate  how  much  the  proposed  material-model  modifica¬ 
tions  can  affect  the  composite  laminate  computed  ballistic  limit, 
as  quantified  by  the  V50  (the  projectile  incident  velocity  at 
which  the  probability  of  target  penetration  is  50%). 

4.1  Uniaxial  Tension  Test 

In  this  section,  a  quasi-static  uniaxial  tension  test  is  simulated 
using  the  conventional  finite  element  method.  Details  of  the 
finite  element  procedure  used  are  identical  to  the  one  reported  in 
our  prior  study  (Ref  11,  20-22).  In  short,  a  rectangular- 
parallelepiped-shaped  composite-material  test  sample  with  its 
longitudinal  axis  aligned  in  the  weft  direction  is  subjected  to  a 
progressively  larger  axial  strain  (while  no  lateral  constraints  are 
imposed,  to  attain  a  uniaxial  stress  state  condition).  To  identify 
differences  in  the  tensile-response  of  the  test  sample  brought 
about  by  the  proposed  modifications  in  the  material  model, 
various  global-response  and  spatial  distributions  of  material 
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Fig.  3  Topological  relationships  between  (a)  fibrils  and  fibers;  and 
(b)  fibers  and  yams  adopted  in  the  present  investigation 


state  variables  are  monitored/generated.  Details  of  the  results 
obtained  are  presented  below. 

A  comparison  of  the  uniaxial  tension  (loading)  true  stress 
versus  uniaxial  tension  true  strain  curves  for  a  prototypical 
Kevlar®  plain  woven  composite  material  obtained  using  the 
original  Yen  (Ref  8)  material  model  and  its  upgrade  proposed  in 
the  present  study  is  shown  in  Fig.  6.  Examination  of  this  figure 
reveals  that 

(a)  almost  up  to  the  onset  of  failure,  the  two  curves  are 
effectively  coincident;  (b)  However,  the  composite  updated 
material  model  shows  that  the  onset  of  failure  occurs  at  a  3-4% 
lower  stress  level  and  a  slightly  smaller  strain;  and  (c)  The 
updated  material  model  predicts  a  higher  rate  of  material 
damage/degradation.  In  other  words,  the  material  displays 
almost  a  brittle-like  behavior. 

In  Fig.  7(a)-(d)  and  8(a)-(d),  a  comparison  is  made  between 
temporal  evolutions  of  spatial  distribution  of  the  axial  tensile 
stress  in  the  weft  direction  in  the  same  prototypical  Kevlar® 
plain  woven  composite  material,  obtained  using  the  original 
Yen  (Ref  8)  material  model  and  the  updated  material  model 
proposed  in  the  present  study.  An  examination  of  these  results 
reveals  that,  at  each  level  of  the  imposed  overall  axial  tensile 
strain,  the  spatial  distribution  of  the  axial  tensile  stress  is  fairly 
uniform.  This  finding  is  somewhat  influenced  by  the  choice  of 
the  stress  scale  (100-600  MPa)  used  in  these  figures.  In  other 
words,  to  obtain  distinct  colors  for  specific  stress  levels  over  the 
entire  loading  cycle,  a  constant/common  stress  range  was  used. 

In  Fig.  9(a)  and  (b),  a  comparison  is  made  between  spatial 
distribution  of  the  axial  tensile  stress  in  the  weft  direction  in  the 
same  prototypical  Kevlar®  plain  woven  composite  material, 
obtained  using  the  original  Yen  (Ref  8)  material  model  and  the 
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Fig.  4  The  effect  of  the  concentration  of  the  (a)  -COOH  chain  end 
group;  and  (b)  -S03H  side  group  on  the  orthotropic  elastic  stiffness 
constants  of  the  Kevlar®  fibrils 


updated  material  model  proposed  in  the  present  study  at  the 
instant  when  the  overall/average  axial  stress  experiences  a 
maximum.  An  examination  of  the  results  displayed  in  these 
figures  shows  that,  while  in  the  case  of  the  original  Yen  (Ref  8) 
material  model,  Fig.  9(a),  spatial  distribution  of  the  axial  stress 
is  indeed  quite  uniform,  in  the  case  of  the  upgraded  composite 
material  model  the  axial-stress  distribution  is  somewhat 
nonuniform.  To  explain  these  differences,  a  comparison  is 
made  in  Fig.  10(a)  and  (b)  between  the  spatial  distributions  of 
the  extent  of  axial  tension/trans verse  shear  damage  in  the  weft 
direction  in  Kevlar®  plain  woven  composite  material  using  the 
original  Yen  (Ref  8)  material  model  and  the  updated  material 
model  proposed  in  the  present  study  at  the  instant  when  the 
overall/average  axial  stress  experiences  a  maximum.  It  is  seen 
that,  while  in  the  case  of  original  Yen  (Ref  8)  material  model 
this  damage  variable  is  fairly  uniformly  distributed  (and  equal 
to  zero),  in  the  case  of  the  updated  material  model,  one  of  the 
ten  cubic  finite  elements  constituting  the  tensile  specimen  has 
acquired  a  non-zero  value  of  this  damage  variable.  This  finding 
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Table  2  A  summary  of  the  effect  of  the  most  common  defects  found  in  Kevlar®-based  materials  and  their  (same  defect- 
type)  clusters  on  the  PPTA  fibril’s  axial  tensile  strength  (in  GPa) 


Flaw  group 

Flaw 

Number  of  defects  in  the  largest  cluster 

1  2  3 

Isolated  chain  ends  (point  defect) 

-COOH 

6.0 

5.6 

5.0 

-COO“Na+ 

-nh2 

6.1 

5.8 

5.2 

-NH3+HS04“ 

Side  groups  (point  defect) 

-so3h 

5.8 

5.1 

4.1 

-SOs“Na+ 

Voids  and  interstitials  (point  defects) 

Microvoids 

3.8 

N/A 

N/A 

Mobile  Trapped  H2S04 

3.9 

N/A 

N/A 

Defect  bands  (planar  defects) 

NH3+HS04  agglomerated  chain  ends 

3.2 

N/A 

N/A 

Table  3  A  summary  of  the  effect  of  the  composition  and 
the  size  of  the  mixed  S03H/S03-Na+based  side  group  and 
COOH/COO_Na+  based  chain  end  defect  clusters  on  the 
PPTA  fibril’s  axial  tensile  strength  (in  GPa) 


Number  of  SO3H/SO3  Na+  defects  in  the 
cluster 


1 

2 

3 

Number  of  COOH/COO“Na+ 

defects  in  the  cluster 

1  5.5 

4.8 

3.8 

2  5.3 

4.4 

3.4 

3  4.5 

3.8 

2.7 

explains  not  only  the  axial-stress  nonuniformities  observed  in 
Fig.  9(b)  but  also  the  more  brittle  nature  of  the  composite 
material.  In  other  words,  once  one  of  the  elements  acquires  a 
non-zero  value  of  the  uniaxial  tension/transverse  shear  damage 
variable,  further  deformation  becomes  localized  within  this 
element.  This  causes  both  the  structural  weakening  of  the 
material  in  this  element  (due  to  loss  of  the  associated  Young’s 
modulus)  and  the  geometrical  weakening  (due  to  a  reduction  in 
the  cross-sectional  area  supporting  the  applied  axial  load).  It 
should  be  noted  that  the  observed  differences  between  the 
results  presented  in  Fig.  9(a)  and  (b),  as  well  as  the  differences 
between  the  results  presented  in  Fig.  10(a)  and  (b),  are  due  to 
the  stochastic  nature  of  material  strength  parameters  imple¬ 
mented  into  the  updated  material  model.  In  other  words,  the 
failed  element  identified  in  Fig.  9(b)  and  10(b)  is  associated 
with  the  lowest  level  of  tensile  strength. 

The  results  displayed  in  Fig.  6-10  clearly  reveal  that 
incorporation  of  the  realistic  micro-mechanics  effects  such  as 
the  operation  of  various  fiber/matrix-damage  mechanisms  and 
the  stochastic  nature  of  the  material  strength  can  significantly 
alter  the  global  as  well  as  the  local  response  of  the  material  to 
the  imposed  loading.  While  this  finding  can  be  considered  as 
expected,  its  direct  quantitative  experimental  validation  is  quite 
challenging  and  is  beyond  the  scope  of  the  present  study. 

4.2  Laminate  Transverse  Impact  by  a  Rigid  Cylindrical 
Projectile 

In  this  section,  the  problem  of  composite  laminate  transverse 
impact  by  a  right  circular  cylindrical  rigid  projectile  is  analyzed 


Fig.  5  (a)  Kevlar®-fibril’s  axial  strength  probability  density  func¬ 

tion;  and  (b)  the  corresponding  cumulative  distribution  function 
determined  in  the  present  study  using  molecular-type  computational 
analyses 


692— Volume  22(3)  March  2013 


Journal  of  Materials  Engineering  and  Performance 


700 


-+ -  Original  Yen  Model 


Fig.  6  A  comparison  of  the  uniaxial  tensile  true  stress  vs.  uniaxial 
tensile  true  strain  behavior  (in  the  weft  direction)  in  a  prototypical 
Kevlar®  plain  woven  composite  material  obtained  using  the  original 
Yen  (Ref  8)  material  model  and  its  upgrade  proposed  in  the  present 
study 
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Fig.  8  Spatial  distribution  of  the  axial  tensile  stress  in  the  weft 
direction  in  a  prototypical  Kevlar®  plain  woven  composite  material 
obtained  using  the  upgraded  model  proposed  in  the  present  study  at 
the  average  axial  strain  levels  of  (a)  0.004;  (b)  0.014;  (c)  0.024;  and 
(d)  0.032 
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Fig.  7  Spatial  distribution  of  the  axial  tensile  stress  in  the  weft 
direction  in  a  prototypical  Kevlar®  plain  woven  composite  material 
obtained  using  the  original  Yen  (Ref  8)  material  model  at  the  aver¬ 
age  axial  strain  levels  of  (a)  0.004;  (b)  0.014;  (c)  0.024;  and  (d) 
0.033 

using  a  dynamic  explicit  finite  element  computational  analysis. 
Since  details  of  the  finite  element  procedure  used  is  identical  to 
the  one  reported  in  our  prior  study  (Ref  11),  they  will  not  be 
repeated  here. 

An  example  of  the  prototypical  results  obtained  in  the 
present  study  is  given  in  Fig.  1  l(a)-(d).  The  results  displayed  in 
Fig.  ll(a)-(d)  pertain  to  the  spatial  distribution  of  the  compos¬ 
ite-laminate  material  during  impact.  It  should  be  noted  that  to 
enable  monitoring  of  the  penetration  process,  the  target  plate  is 
made  transparent.  The  results  displayed  in  Fig.  1  l(a)-(d)  clearly 
reveal  the  development  of  the  back  face  bulge.  In  general,  it 
was  found  that  details  regarding  the  spatial  distribution  and  the 
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Fig.  9  A  comparison  between  the  spatial  distributions  of  the  axial 
tensile  stress  in  the  weft  direction  in  Kevlar®  plain  woven  composite 
material  at  the  instant  when  the  overall/average  axial  stress  experi¬ 
ences  a  maximum,  obtained  using  (a)  the  original  Yen  (Ref  8)  mate¬ 
rial  model;  and  (b)  the  updated  material  model  proposed  in  the 
present  study 
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Fig.  10  A  comparison  between  the  spatial  distributions  of  the  ex¬ 
tent  of  axial  tension/transverse  shear  damage  in  the  weft  direction  in 
Kevlar®  plain  woven  composite  material  at  the  instant  when  the 
overall/average  axial  stress  experiences  a  maximum,  obtained  using 
(a)  the  original  Yen  (Ref  8)  material  model;  and  (b)  the  updated 
material  model  proposed  in  the  present  study 

temporal  evolution  of  the  target  material  during  impact  are 
affected  by  the  modifications  made  to  the  original  Yen  (Ref  8) 
composite  material  model.  Details  regarding  the  effect  of  these 
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Fig.  11  Typical  results  pertaining  to  the  temporal  evolution  of  the 
spatial  distribution  of  the  composite-laminate  material  during  impact 
by  a  right-circular  cylindrical  rigid  projectile 
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Fig.  12  A  typical  projectile  residual  velocity  vs.  projectile  incident 
velocity  plot  showing  differences  in  V50  as  predicted  by  the  original 
Yen  (Ref  8)  material  model  and  its  upgrade  proposed  in  the  present 
study 


material  model  modifications  on  the  temporal  evolution  and 
spatial  distribution  of  the  target-material  and  its  various 
mechanical  and  microstructural  fields  will  be  presented  in  our 
future  communication.  In  the  remainder  of  this  section,  results 
pertaining  to  the  effect  of  material  model  modifications  on  V50 
are  presented  and  discussed. 


Typical  results  pertaining  to  the  effect  of  projectile  initial 
velocity  on  the  projectile  residual  velocity  in  the  case  of  the 
original  Yen  (Ref  8)  composite  material  model  and  its  present 
modification  are  displayed  in  Fig.  12.  It  is  seen  that  V50 
(defined  here  as  the  projectile  initial  velocity  at  which  the 
projectile  residual  velocity  is  zero)  predicted  by  the  analysis  in 
which  the  modified  composite  material  model  is  used  is  ca.  5% 
lower  than  its  counterpart  associated  with  the  use  of  the  original 
Yen  (Ref  8)  composite  material  model.  This  effect  is  deemed  as 
being  significant  and  justifies  the  present  efforts  aimed  at 
incorporating  fine-scale  microstructural  effects  into  the  contin¬ 
uum-level  material  model.  However,  no  relevant  experimental 
data  are  presently  available  to  be  able  to  judge  the  potential 
improvements  in  the  experiment/calculation  agreement  brought 
about  by  the  use  of  the  modified  composite  material  model. 
This  deficiency  in  the  relevant  experimental  data  will  be 
addressed  in  our  future  study. 


5.  Summary  and  Conclusions 

Based  on  the  results  obtained  in  the  present  study,  the 
following  summary  remarks  and  main  conclusions  can  be 
drawn: 

1 .  A  continuum-level  rate-dependent  material  model  for 
plain  woven  Kevlar®-fiber-reinforced  polymer  matrix 
composite  has  been  upgraded  by  incorporating  the  effects 
of  fibril/fiber/yam  level  microstructure  phenomena  and 
processes. 

2.  Specific  microstructural  effects  incorporated  into  the  con¬ 
tinuum-level  material  model  account  for  the  role  of  vari¬ 
ous  material  synthesis/fiber  processing-induced  defects  on 
the  material  strength,  or  more  specifically  on  the  material 
strength  distribution  function. 

3.  For  a  prototypical  plain  woven  Kevlar®-fiber-reinforced 
polymer  matrix  composite,  the  results  obtained  using  a 
series  of  finite  element  analyses  show  that  the  predicted 
response  of  the  material  both  under  quasi-static  and 
dynamic  ballistic-impact  conditions  is  significantly  affected 
as  a  result  of  the  incorporation  of  microstructural  effects. 
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